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Abstract: In wireless sensor networks, data aggregation routing could reduce the number 
of data transmissions so as to achieve energy efficient transmission. However, data 
aggregation introduces data retransmission that is caused by co-channel interference from 
neighboring sensor nodes. This kind of co-channel interference could result in extra energy 
consumption and significant latency from retransmission. This will jeopardize the benefits 
of data aggregation. One possible solution to circumvent data retransmission caused by co- 
channel interference is to assign different channels to every sensor node that is within each 
other's interference range on the data aggregation tree. By associating each radio with a 
different channel, a sensor node could receive data from all the children nodes on the data 
aggregation tree simultaneously. This could reduce the latency from the data source nodes 
back to the sink so as to meet the user's delay QoS. Since the number of radios on each 
sensor node and the number of non-overlapping channels are all limited resources in 
wireless sensor networks, a challenging question here is to minimize the total transmission 
cost under limited number of non-overlapping channels in multi-radio wireless sensor 
networks. This channel constrained data aggregation routing problem in multi-radio 
wireless sensor networks is an NP-hard problem. I first model this problem as a mixed 
integer and linear programming problem where the objective is to minimize the total 
transmission subject to the data aggregation routing, channel and radio resources 
constraints. The solution approach is based on the Lagrangean relaxation technique to relax 
some constraints into the objective function and then to derive a set of independent 
subproblems. By optimally solving these subproblems, it can not only calculate the lower 
bound of the original primal problem but also provide useful information to get the primal 
feasible solutions. By incorporating these Lagrangean multipliers as the link arc weight, 
the optimization-based heuristics are proposed to get energy-efficient data aggregation tree 
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with better resource (channel and radio) utilization. From the computational experiments, 
the proposed optimization-based approach is superior to existing heuristics under all tested 
cases. 

Keywords: channel assignment; data aggregation; multi-radio; optimization; wireless 
sensor networks 



1. Introduction 

With the capability of sensing, computing and communication embedded on the sensor node, the 
wireless sensor network (WSN) is a promising technology to probe and collect environmental 
information. Without the necessity of expensive wiring cost for constructing the sensor network, the 
WSN could deploy sensor nodes at any location more efficiently [1]. A representative WSN is shown 
in Figure 1. Sensor nodes are usually scattered in a sensor field. When any event occurs, such as 
surging irradiation or temperature declining below certain threshold, sensor nodes within specific 
sensing range (data source nodes) detect this event and collect data which would be transmitted to the 
sink node for taking further processing. 



The application scenario described above is called event-driven that sensors are assigned to detect a 
particular event, which is shown in Figure 1(a). There are two other different applications of wireless 
sensor networks, which are namely periodic and query-based. These two applications could be 
classified as random-source, which is shown in Figure 1(b). Periodic scenario sensors probe 
environmental information periodically and report their measurements back to the sink node. All 
sensors in this kind of networks are necessitated to be synchronized such that all sensors sense 
information and report it simultaneously. Query-based scenario is applied to user-oriented applications. 
User can query information from certain area of sensors to require measurements that he is interested 



Figure 1. Typical wireless sensor networks. 




(a) Event driven model 



(b) Random source model 
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in. As can be seen in Figure 1(b), the data source nodes in random source model are not clustered as in 
event driven model. 

Power efficient communication in the WSN is an interesting and blooming research area [2]. It is 
almost impossible to replace the battery in the sensor node due to the limited energy power in the 
sensor node. In traditional wireless network, the data transmission is unicasting (one to one) or 
multicasting (one to many). In the WSN, multiples data source nodes sense the data and send them 
back to only one sink node. Hence, the data communication in the WSN is a kind of reverse 
multicasting {many to one), also known as data aggregation routing. This makes power efficient 
communication in the WSN different from traditional wireless network. 

For data aggregation routing, raw data from multiple children sensor nodes are collected and 
processed before transmission. Data aggregation can minimize the number of transmission by 
eliminating redundant data from different source nodes [3-5]. Figure 2 gives an example of data 
aggregation routing and address centric routing where the maximum temperature is reported to a sink 
node. Label x(y) at each node represents the local temperature measurement is x while the aggregated 
(maximum) value so far is y. For example, at node 27(30), the maximum temperature up to now is 
MAX (27, 30) = 30. Assume the transmission cost on each link is equal to 1.0. In traditional address 
centric routing shown in Figure 2(b), each Origin-Destination (OD) pair follows the shortest path. 
Then the total cost for these three OD pairs is 4.0. However, with the capability of data aggregation at 
each sensor node, a more power efficient transmission is shown in Figure 2(a) via data aggregation 
routing where the total cost is 3.0. 

Figure 2. Data aggregation in data aggregation routing and address centric routing. 
30(30) 28 30 28 




(a) Data aggregation routing (b) Address centric routing 

Although data aggregation in the WSN could reduce the number of transmissions to save 
transmission cost, it could introduce additional MAC layer retransmission energy loss. Based on the 
CSMA/CA protocol, data transmission from multiple sensor nodes to the same sensor node for data 
aggregation will incur collision. When there is collision, retransmission is required to ensure that the 
data is successfully received at the aggregation node and this will incur additional energy consumption. 
Basically, the more flows are aggregated at the sensor node, the higher probability that the senders will 
incur data retransmission. In Figure 3(a), a node ns (with three children nodes) will suffer severe 
collisions which results in more retransmission times as compared to a node ns (with two children 
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nodes) in Figure 3(b). Besides extra energy consumption, retransmission also incurs additional latency, 
which is unacceptable in delay sensitive applications. The extra energy consumption and additional 
latency from retransmission will jeopardize the advantage from data aggregation. 

Figure 3. Channel, radio and rata aggregation routing. 
Data source Node: n ls n 2 , n 3 




Total Channels Required = 4 Total Channels Required = 3 

Number of Radios Required = 3 Number of Radios Required = 2 

Total Transmission Cost = 21.84 Total Transmission Cost = 25.68 

(a) Pure data aggregation routing (h) Channel and Radio Aware 



data aggregation routing 

By assigning different channels to the sensor nodes that are within each other's interference range, 
the retransmission problem caused by collision could be circumvented. If there is sufficient number of 
channels, then we could assign a different channel to every sensor node on the aggregation tree such 
that there is no extra energy loss from retransmission. In the meantime, the latency could also be 
minimized. However, the number of channels is a limited and valuable resource in wireless networks. 
For instance, in IEEE 802.11b, there are only three non-overlapping channels [6]. Hence, the question 
is how to assign the limited channels to the sensor nodes on the aggregation tree such that the total 
transmission power could be minimized. 

Besides limited number of available non-overlapping channels, number of radios on each sensor 
node is also a limited resource. If two children sensor nodes use two different channels transmit data 
back to the same sensor node, then this sensor node will need two radios to receive data 
simultaneously. Otherwise, it will incur larger latency for a single radio sensor node to switch different 
channel to receive from its children nodes. Hence, from the latency point of view, for any sensor node 
that is on the data aggregation tree, the number of radios equipped on this node must be greater than or 
equal to the number of children nodes. In Figure 3, 1 illustrate an example where the sensor nodes are 
randomly placed in a 15 x 15 area. In this example, the transmission cost is equal to the square of the 
Euclidean distance of the transmission radius. The sensor nodes with the same color are assigned with 
the same channel (e.g., ni and n^). To prevent collision, any two sensor nodes that are within each 
other's transmission range could not assign the same channel (e.g., «2 and «j). In addition, even though 
«3 and n4 are not within each other's transmission range, n 3 still could not reuse n4 's channel because of 
«j. This is referred to as the hidden node problem. In this case, we say n 3 and n^are within each other's 
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interference range. It is important to note that interference range is larger than the transmission range 
because of the hidden node problem. If every sensor node has the same transmission radius, then 
interference range is equal to the twice of transmission range in order to capture the co-channel 
interference from hidden nodes. In Figure 3(a), when the number of channels and radio is unlimited, 
the minimum total transmission cost is 21.84. The total number of channels required is 4 and the 
number of radios required for a node n 5 is 3. If the total number of channels is limited to 3 or the 
number of radios on each node is limited to 2, then Figure 3(a) is not a feasible solution. In Figure 3(b), 
it shows the data aggregation tree for channel and radio aware data aggregation routing when the total 
number of channels is limited to be 3 and the number of radios on each sensor node is limited to 2. The 
transmission cost is 25.68, which is slightly larger than Figure 3(a). 

To perform Channel and Radio Constrained Data Aggregation Routing (CRDAR) in the WSN is 
even more challenging than pure data aggregation routing in the WSN. The channel assignment in 
wireless network could be modeled as a graph coloring problem in graph theory where adjacent nodes 
could not be assigned with the same color. This graph coloring problem is proven to be a NP-hard 
problem [7]. CRDAR that contains the channel assignment problem is also an NP-hard problem. In 
this paper, for the first time, I first model the CRDAR problem as a mixed integer linear programming 
(MILP) problem. Then I relax some of the constraints and derive a set of independent subproblems. By 
optimally solving each independent subproblem and adjusting the Lagrangean multipliers at each 
iteration by the subgradient method, we get the tightest lower bound of the CRDAR problem. By 
utilizing the solution to the Lagrangean dual problem and the information from the Lagrangean 
multipliers, a getting primal heuristic algorithm (LGR) is devised to identify the channel and radio 
constrained minimum transmission cost data aggregation tree. 

Note that this integrated channel assignment and routing problem is also an important issue in 
multi-radio wireless networks. Hence, the proposed optimization-based algorithm could also be 
applied to general wireless networks. The reason that I specifically focus on the WSN is because the 
nature of "many-to-one" communication from multiple data source nodes to one sink in the WSN is 
different from "one -to-many" communication in the wireless networks. This many-to-one 
communication will increase the probability of co-channel interference so as to make the integrated 
channel assignment and routing problem in the multi-radio WSN more challenging than the traditional 
wireless network. 

The remainder of this paper is organized as follows. In Section 2, existing literature on data 
aggregation routing and channel assignment problem in wireless networks is surveyed. In Section 3, I 
formulate the CRDAR problem as the MILP mathematical problem. In Section 4, Lagrangean 
relaxation scheme is applied to relax some constraints and algorithms are proposed to solve the 
Lagrangean dual problem optimally. In Section 5, the novel optimization-based heuristics are devised 
to get the primal feasible solution. In Section 6, the numerical results and performance comparisons 
are demonstrated. Finally, concluding remarks are summarized in Section 7. 

2. Related Works 

In wireless network, if the transmission radius of a node is r, then the power consumption is 
measured as r a + c, where a is a signal attenuation constant (usually between 2 to 4) and c is a positive 
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constant that represents signal processing [3]. In [8], they study the tradeoff between power 
consumption transmission radius and coverage of the transmission node. For long transmission radius, 
more sensor nodes could be covered so that the total number of transmission could be reduced. 
However, long transmission radius incurs significant power consumption, especially for large a, that 
would sacrifice the gain from reduced total number of transmission. 

Pure data aggregation routing problem in the WSN has been studied by several works. In [8], they 
propose centralized heuristic based on Prim's shortest path algorithm to construct a data aggregation 
tree. However, as shown in [9], the data aggregation tree constructed by shortest path algorithm 
(shortest path tree, SPT) does not facilitate the data aggregation advantage. In [5], three interesting 
suboptimal aggregation heuristics, Shortest Paths Tree (SPT), Center at Nearest Source (CNS), and 
Greedy Incremental Tree (GIT), are proposed. It is shown that GIT could get the best results. The idea 
of GIT scheme is initially the only member in the tree is the sink node. Each data source find the 
shortest hop path to this tree and the data source with the minimum hop along with the intermediate 
nodes on this path are included in this tree. This process is repeated until all source nodes are included 
in the tree. 

In [3,9], they propose a rigorous MILP formulation for data aggregation routing problem and 
propose solution approaches based on Lagrangean relaxation. Optimization-based heuristic are 
proposed to solve the pure data aggregation routing problem. From the computational experiments, 
their optimization-based algorithm is superior to the SPT, CNS and GIT [5] in both random-source 
model and event-driven model. In [4], they consider the MAC aware data aggregation routing in the 
WSN. They capture the energy consumption tradeoffs between the data aggregation and retransmission 
in the CSMA/CA MAC protocol by proposing an interesting optimization-based heuristics. It is shown 
that proposed algorithms could construct more energy efficient data aggregation tree with considering 
MAC layer retransmission mechanism than existing data aggregation algorithms. 

In [11], the channel assignment problem is modeled as a graph coloring problem and proposes an 
optimization-based heuristics to tackle this issue. However, the hidden node problem is not addressed 
in [11]. Joint channel assignment and routing in multi-hop wireless networks has been studied 
in [12,13]. They address the interference problem in wireless network by jointly channel assignment 
and routing algorithm in order to achieve maximum throughput/system capacity. However, power 
efficient communication is not addressed in [12,13]. In addition, they addressed the unicasting routing 
problem which is not applicable to data aggregation routing problem in the WSN. 

In [14], they study the tradeoff between data aggregation and latency in the WSN. Data aggregation 
tree is constructed by using the earliest-first, randomized, nearest-first and weighted-randomized to 
identify the parent node to relay the data from the data source node back to the sink node. Basically, by 
assigning different time slot (i.e., channel) to every sensor node on data aggregation tree that has the 
same parent node, there will be no collision but introducing large latency. Hence, there is a tradeoff 
between data aggregation and latency. In [15], they consider the latency issue in constructing a 
minimum energy data aggregation tree. A data aggregation tree is a balanced binary tree where 
initially the sink node finds the nearest two sensor nodes as its children, and each children node 
identifies another two nearest nodes as its children node. This process is repeated until all data source 
nodes are included in this balanced data aggregation tree. After the data aggregation tree is determined, 
channel assignment is performed to minimize latency and transmission power. This two-phase 
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approach is also shown in [16]. However, restricting the data aggregation tree to be the balanced 
binary tree might lead to long data aggregation tree that has larger transmission cost. 

In [17], they show how to minimize energy consumption by using the prioritized channel 
assignment and scheduling the listen and sleep time. Sensor nodes with less residual power are 
assigned with higher priority for channel assignment. Sensor nodes with low priority are scheduled to 
sleep to save energy loss. However, data aggregation routing is not addressed in [17]. In [18], sensor 
networks are partitioned into clusters and each sensor node sends its data to its cluster head before 
transmitting outside the cluster. All inter-cluster communication is done by the cluster head. However, 
channel reuse for better channel utilization is not addressed in [18]. 

In [19], they propose four interesting heuristic algorithms (SPT, GIT, CAGIT and ICADAR) to 
solve the channel constrained data aggregation routing problem. The first three heuristics are first to 
identify the data aggregation tree and then perform channel assignment to satisfy the channel 
constraint, which is a one-shot algorithm. The fourth algorithm is is an iterative procedure where 
routing and channel assignment may perform several times to identify a feasible and energy efficient 
data aggregation tree. However, the radio constraint is not addressed in [19] such that it might not be 
feasible in the radio-constrained WSN. 

3. Problem Formulation 

The CRDAR is formulated as a MILP problem. The objective function is to minimize the total 
transmission cost. The constraints include the data aggregation tree, co-channel interference constraint, 
and channel and radio resource constraint. I consider the case where multiple events occur 
simultaneously and send back to the sink node via different data aggregation trees and each event is 
modeled as one multicast group. Each event carries different data such that data aggregation could not 
be performed between different data aggregation trees. If any sensor node is on two different data 
aggregation tree, then this sensor node needs two channels/radios to transmit the data on different tree 
respectively. Hence, more channels/radios will be needed as compared to single data aggregation tree. 

3.1. Notations in the Formulation 

The notations used in the formulation are as follows. 
Input variables: 

N : the set of sensor nodes in the WSN; 

L : the set of possible communication links in the WSN; 

G : the set of multicast groups; 

D g : the set of data source nodes of multicast group g e G; 

F : the set of non-overlapping channels in the WSN; 

ai : unit power transmission cost on link / e L; 

P g d : the set of candidate paths from the sink node of multicast group g to its data source node d; 

h g : The minimum number of hops to the farthest destination node in multicast group g; 

S p i : = 1, if link / is on the path p; = 0, otherwise; 
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Sjk : = 1, if a sensor node j can not use the same channel with a sensor node k; = 0, otherwise; 
pik : = 1, if a sensor node k is the termination node of link /; = 0, otherwise; 
oik : = 1, if a sensor node k is the source node of link /; = 0, otherwise; 
co : the number of available non-overlapping channels, i.e., co = \F\; 
Rj : the number of radios equipped on a node j e N. 
Decision variables: 

x g pd : = 1, if the sink node of multicast group g use path p to reach its data source node d; = 0, 
otherwise; 

= 1, if link / is on the data aggregation tree of multicast group g; = 0, otherwise; 
aggregate flow on link /; 

= 1, if channel i is assigned to a data source node j; = 0, otherwise; 
= 1, if channel i is assigned to any data source node; = 0, otherwise. 
Basically, the size of P g d will grow exponentially with number of sensor nodes (i.e., |iV|.). So it is 
almost impossible to enumerate all possible paths at large network size. I will show in Sec. 4 that we 
do not need to enumerate all possible paths for P gc i. The Lagrangean multipliers associated with 
decision variable x gp d enable us to identify the shortest path for every data source node d of multicast 
group g by using the Dijkstra's shortest path algorithm. Hence, unlike commercial optimization 
package (e.g., CPLEX) that we need identify all the values for P gd , P gd is just a notation for our 
proposed Lagrangean relaxation scheme. 

Another important parameter is ey*, which captures the co-channel interference. Recall from the 
example shown in Figure 3, interference range is larger than the transmission range due to the hidden 
node problem. Without loss of generality, I assume wireless links are symmetrical and the interference 
range is twice of the transmission range. In other words, if a sensor node k is a two-hop neighboring 
node for a sensor node j, then Sjk = 1 and % = 1 . 



3.2 MILPfor CRDAR Problem 



Problem (P): 

Objective function: 



subject to 
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LA XgpdSpl - VgeG,/eZ (6) 



S Xgpd = i Vg e G , J e D 



(7) 

- Oorl VpeP gd , geG, d^D g (g) 

YuYy s iPij < H m y VjeN (9) 

SZV// < Z«, 7 V/eJV (10) 

(my + m ik )e Jk < { VieF,jeN,keN (11) 

my = Oorl VieFJeN (12) 

< », Vie F ,j g N (13) 
n, = Oorl VieF (14) 

< co (15) 

S w «7 < Rj V/eiV. (16) 

The objective function is to minimize the total transmission cost. For instance, in Figure 4, there are 
two data aggregation trees for two events and each event is referred to as one multicast group. So at a 
sensor node n?, it needs to send two data (one from multicast group 1 and the other from multicast 
group 2) back to the sink node. Then on link / between a node ri7 and a node its, Q=2. However, on 
link / between a node ns and a node nj, Q = 1, because the data are aggregated for the same multicast 
group at a node n s . The unit power transmission cost on link (a/) here is if, where d is the Euclidean 
distance of link / and a is the signal attenuation constant (usually between 2 to 4). By this definition, 
the total cost of the reverse multicast tree from the data source nodes back to the sink node is identical 
to the total cost of the multicast tree from the sink node to the data source nodes. Hence, the total 
transmission cost for multicast tree (e.g., Figure 4(b)) is identical to the data aggregation tree (e.g., 
Figure 4(a)) except the transmission is in opposite direction. By this observation, in the following 
constraints, I constructs a minimum cost multicast tree for each multicast group from the sink node to 
all its data source nodes and in the meantime to satisfy the co-channel interference, channel and radio 
resource constraints. 

Note that in considering the total power consumption, energy consumption in the idle mode is 
significant such that the sleep/awake mechanism for sensor nodes plays an important role to minimize 
the total power consumption. In this paper, I only addresses the transmission cost instead of total 
power consumption for the CRDAR problem. Therefore, the sleep/awake mechanism is outside the 
scope of this paper. 
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Figure 4. Data aggregation trees and the corresponding multicast trees in the WSN. 




Data source nodes for 
multicast group 1: ni, n 2 , n 3 

Data source nodes for 
multicast group 2: n 4 , n 5 , n 6 
selected link on 
^ group 1 

selected link on 
group 2 

' Q 

^ . Event source l - 
Sink Node: S 



A' Sensing 
range 




(a) data aggregation tree 



(b) multicast tree 



Constraints (1), (2) and (4) enforce that the aggregate flows on link / must be larger than or equal to 
the number of data aggregation tree adopts link /. In Figure 4(b), link rig—My is included in two 
multicast trees, then C„ s ^„ 7 must be at least two. Since the objective function is to minimize the Q, 

then C n will be equal to two at the optimal solution. Hence, it should be an equality at 

Constraint (1), by changing it to smaller than or equal to is a kind of relaxation. From the above 
argument, the equality will be hold at the optimal solution. Hence, the data aggregation property in the 
WSN is also implicitly enforced in these three constraints. 

In Constraint (3), it enforce that if any multicast group g adopts a path to reach its data source node 
d, then this path must be included by the multicast tree of multicast group g. Constraint (5) enforces 
that the number of links on the data aggregation tree of multicast group g must be at least h g and \D g \. 
Hence, h g and D are the legitimate lower bound. For example, in Figure 4, there are three data 

source nodes for multicast group 1, then D =3. The minimum hop count for each data source node 



could be obtained by running the Bellman-ford algorithm. Then minimum hop count for m, ri2 and ni 
is 4, 5 and 3 respectively. In this case, h g = 5. Then the number of selected links (i.e., T,y g} ) on the 

multicast tree of multicast group 1 must be at least 5. 

On the left hand side of Constraint (6), it calculates the number of paths that traverse over link / for 
each multicast group. This number is at most \D g \. If there is a cycle on the union of the routing paths 
to every data source node and link / is on the cycle, then Constraint (6) will not be satisfied because 
there will be unlimited number of paths traversing over link /. Hence, Constraint (6) is to enforce that 
there will be no cycle on the union of the routing paths to every data source node. Constraints (7) 
and (8) enforce that the sink node selects exactly one routing path to every data source node. In 
addition, the objection function is to minimize the total transmission cost, so more than one incoming 
link traversing over a sensor node will incur larger cost. This guarantees that every sensor node has at 
most one incoming link. Hence, Constraints (6), (7) and (8) in conjunction with the objective function 
enforce the union of the routing path to every data source node of the multicast group shall be a tree. 
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For example, in Figure 4, the multicast tree for multicast group 1 is the union of routing paths to three 
data source nodes (i.e. S—>n9—>n?j—>n7—>ni, S—>n9—>n8—>n7—>ni—>n2, S—>n9—>ti8—>n7—>n3). 

In Constraint (9), it enforces that the number of assigned channels on a sensor node must be no less 
than the number of multicast groups selecting this node to transmit data. Recall that the data 
aggregation tree is the same as multicast tree except in opposite transmission direction. The number of 
multicast links terminated at this node indicates the number of multicast groups include this node on 
the routing path to transmit data back to the sink. For example, in Figure 4(b), a sensor node ri2 is the 
termination node of link ni—>ri2, so a node ^must be assigned with a channel to transmit data back to 
fii. Sensor node ri7 is the termination node of link ns—> n 7 , and link n$—> n 7 is selected by both group 1 
and 2, so a node n 7 must be assigned with two channels to transmit data to each group respectively. In 
Constraint (10), it enforces that the number of assigned channels on a sensor node must be no less than 
the number of its children nodes. This is to insure that the sensor node is assigned with enough 
channels to receive data from its children nodes. For example, in Figure 4(b), a sensor node n 7 is the 
source node of three selected links (n 7 ^n s , n 7 ^nj, n 7 ^n 3 ). Because of the hidden node interference, 
these three nodes (nj, n 3 and ns) will be assigned with different channel to transmit data back to a node 
n 7 . So a node nzmust also be assigned with three channels to receive. In summary, Constraint (9) 
enforces the channel assignment for transmitting and Constraint (10) enforces the channel assignment 
for receiving. 

Constraint (11) is the co-channel interference constraint. When e# = 1 (i.e., nodes j and k are within 
each other's interference range), we can not assign the same channel (say channel i) to nodes j and k. 
In this case, my and mn can not be equal to 1 at the same time. Recall that interference range affects 
two-hop neighbors. Hence, for a sensor node j, every within two hops neighboring node k has ey* = 1 . 
In this case, nodes j and node k could not be assigned with the same channel. For example, in Figure 4, 
a node ri4 is a two-hop neighbor of n& so they can not be assigned with the same channel because of nj. 
Hence, the hidden node co-channel interference is also addressed in this constraint. In Constraint (15), 
it enforces that the total number of used channels should not exceed a>. In Constraint (16), it enforces 
that the number of channels assignment on any node should not exceed the number of radios equipped 
on this node. Hence, radio resource constraint is enforced in Constraint (16). 



In problem (P), there is total number of (\G\ + l]L\ + (\N\ + + £ Z \ p gi 



decision variables. 



Besides extremely large number of decision variables, problem (P) contains channel, radio assignment 
and multicast routing problem. Since channel assignment problem is a coloring problem in the graph 
theory, which is proven to be a NP-hard problem, this makes problem (P) also be a NP-hard problem. I 
propose an optimization-based heuristic based on Lagrangean relaxation scheme to tackle this problem. 



4. Solution Approach — Lagrangean Relaxation 



Constraints (1), (3), (6), (9), (10), (11) and (13) in (P) are relaxed to get the following Lagrangean 
relaxation problem, where ( u) ~ ul ) represents the Lagrangean multiplier. Basically, the more 

constraints are relaxed, the looser duality gap between the solutions to the dual problem and the primal 
problem. Loose duality gap might indicate that the solution to the primal problem might be too far 
from the optimal solution. On the other hand, if too little constraints are relaxed, we might not be able 
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to solve the Lagrangean dual problem optimally. Then the solution to the dual problem is not the true 
lower bound of the primal problem. As I will show in the following paragraph that by relaxing these 
seven constraints in (P) guarantee that the dual problem is optimally solved and in the meantime to 
obtain a tighter duality gap. 



Problem (LR): 



Z LR = min Yfli c i + 2>/ 2>*/-C,| + Z Z E M i 



leL 
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geGdeD leL 
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subject to Constraints (2), (4), (5), (7), (8), (12), (14), (15) and (16). (LR) is decomposed into the 
following five independent subproblems. 
Subproblem 1 : for C, 

minZ(«/-"/)Q (SUB1) 

leL 

subject to (2). 
Subproblem 2: for_y g /. 

min X Tu)y gl - Z Z 5»„ - ZZ^ g |^ + ZZZ"> g ,Ay+ ZZZ<W (SUB2) 

Zel geG geGdeD leL geGleL geGleL jeN geGleL jeN 



subject to (4) and (5). 
Subproblem 3: forx g ^. 

min Z Z Z Z fc/ 

geGdeD g leL peP gd 



pl 



(SUB3) 



subject to (7) and (8). 
Subproblem 4: for m,y. 



min- Z Z("/ + " y 5 K + Z Z Z "M m «7 + m ik + Z Z 



jeNieF 



ieF JeNkeN 



ieF jeN 



(SUB4) 



subject to (12) and (16). 
Subproblem 5: for rij. 



min - X H u l n i 

ieF jeN 

subject to (14) and (15). 



(SUB5) 
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Subproblem 1 is to determine decision variable Q. It can be further decomposed into \L\ 
independent problems. For each link / e L, when the coefficient of C/ (i.e. a r -u) ) is positive, let 

Ci = 0. When the coefficient of Q (i.e. a l -u) ) is negative, let Q= \G\. The computational complexity 

for this algorithm is 0(1) for each link. 

Sub-problem 2 is to determine decision variable y g /. It can be further decomposed into \G\ 
independent problems. For each multicast group g e G, 

Zu]y gl -zf«i \D g |+ 1 u] dl k+ZZ(«;^+«>>„ (SUB2.1) 



mm " 

/el " /el 



/el yeAT 



subject to: 

^ ; =0 or 1 V/gI and J> , ^ max{^, \D g \ } 

/ei 

The algorithm to optimally solve (SUB2.1) is shown as follows. 

Stepl. For every link / where its termination node is sensor node j, the coefficient for y g i on this link / 
is u) - u gl | D | - ^ u 2 gdl + u A j . For every link / where its source node is the sensor node j, the 

d*D g 

coefficient for y g i on this link / is u) - u gl \ D g \ - £ u 2 gdl +u 5 j .Then calculate the number of links whose 
coefficient is negative. 

Step2. For each multicast group ge G, if the number of negative coefficient links is greater than ore 
equal to max{h g , \D g \}, then for these negative coefficient links, assign the corresponding y g i= 1 and let 
other y gi = 0. 

Step3. If the number of negative coefficient links (assume the number is 6) is smaller than 
max{/z g , \D g \), assign the corresponding y g i = 1 for these negative coefficient links. Then sort those 
links that have positive coefficient in ascending order. Assign {max{h g , \D g \}-6} number of smallest 
positive coefficient and let y g i = 1. Finally, let the other y g i = 0. 

The computational complexity of above algorithm is 0(|Z(|Z) g |+log|Z|)) for each multicast group. 

Subproblem 3 is to determine decision variable x gpc i. It can be further decomposed into X|^ g | 

independent shortest path problems with non-negative arc weights {u gdl + u gl ) . They can be effectively 

solved by the Dijkstra's algorithm. The computational complexity of the Dijkstra's algorithm is 0(|A^| 2 ) 
for each data source node of the multicast group. 

Sub-problem 4 is to determine decision variable my. It could be further decomposed into \N\ 
independent subproblems. For each node j e N, 

min - +u 5 j )m ij + £ £ («£ + u% )e jk m y + £ ufa (SUB4.1) 

ieF i^F keN ieF 

subject to: 

my = 0 or 1 V i e F (16) 



Sensors 2009, 9 



4779 



(17) 



In the objective function of (SUB4.1), the coefficient for my is ul - u* - u 5 j + 




Because of Constraint (17), at most Rj number of channels could be selected for a node j. If the number 
of negative coefficient is smaller than Rj, then assign my = 1 for those channels that have negative 
coefficient. If the number of negative coefficient is greater than or equal to Rj, then assign my = 1 for 
those Rj channels that have smallest negative coefficient. Then assign my = 0 for the rest of the 
channels. The computational complexity is 0(\F\\og\F) for each node j e N. 

Sub-problem 5 is to determine decision variable For each channel i e F, since the Lagrangean 
multiplier Uy is positive, so the coefficient of n, (i.e. - Y* u t ) i s negative. Constraint (14) enforce 



that we can at most select co number of m such that n, • = 1. First, we initialize every m = 0. Then by 
sorting these coefficient for each in ascending order, we select co number of n t whose correspond 
coefficient are the smallest and let these n, ■ = 1. In this way, we could have minimum objective value in 
(SUB5). The computational complexity is 0(\F\log\F). 

Applying the above algorithms, we can solve the Lagrangean dual problem (LR) optimally. 
According to weak duality theorem, that is, the objective value of LR is a legitimate lower bound to 
the original problem (P). One can calculate the tightest lower bound and solve the dual problem by 
using the subgradient method [20]. Note that the solutions to the dual problem may not be feasible to 
the primal problem due to several constraints are relaxed. In the sequel, I propose the heuristic for 
getting the primal feasible solution. 

5. Obtaining Primal Feasible Solutions 

The basic idea of getting primal feasible solution (LGR-Primal) is first to identify the energy 
efficient data aggregation and then adjust the routing path to meet the channel/radio resource 
constraint. This is a kind of iterative algorithm which is particularly useful under stringent resource 
constraint (i.e., limited channels and radios). In order to facilitate this idea, the LGR-Primal algorithm 
is first to identify efficient data aggregation tree by using the GIT algorithm. Then perform the channel 
and radio assignment algorithm. If the channel and radio constraint is satisfied, then report the data 
aggregation tree. Otherwise, identify another data aggregation routing path such that the channel and 
radio constraint could be satisfied, this process is repeated until feasible data aggregation tree is 
identified. 

In Figure 5, the data aggregation tree is first determined by the GIT algorithm, and then performs 
channel/radio assignment for every sensor node on the data aggregation tree. The channel/radio 
assignment order has significant impact on the solution quality. According to [19], the data source 
node with shorter hop count on the data aggregation tree to the sink is more likely to aggregate data 
than the data source node with longer hop count. Then by first assigning the channel to the data source 
node with shorter hop count, the data source node with longer hop count could reuse the channel more 
efficiently without violating the co-channel interference constraint. Leverage on this observation, at 
Step 3 of Figure 5, the data source node with shorter hop count to the sink has higher priority for 
channel and radio assignment than the data source node with longer hop count. 
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Figure 5. LGR-Primal Algorithm. 



1 . Assign the arc weight to every link on WSN 

i 

2. Perform GIT algorithm to the identify shortest path from the sink 
node to the every data source node 



T 



3. Assign collision-free channel to every data source node 



F: False 
T: True 




5. identify the sensor node that violate the channel/radio constraint 
and assign the arc weight on every link incident to this node to be Z 



ary. 
oft 



6. Identity the new routing path of the data source node that traverse 
this sensor node by the GIT algorithm 




8. Could not identify feasible solution 



9. Report feasible solution 



At Step 5, it specifies that if we could not assign feasible channel or radio to any sensor node along 
the routing path of the data source node, then we identify another routing path which bypasses the 
sensor nodes that violate the channel or radio constraint for this data source node. By assigning the 
very large arc weight (Z) to the links incident to these sensor nodes, the data source node could bypass 
these sensor nodes by running GIT algorithm again. If the total arc weight for the data source node 
exceeds Z, then we could conclude that there is no feasible channel/radio constrained routing path for 
this data source node. Note that, by assigning the very large arc weight (Z) to the links incident to the 
sensor nodes that violate the channel/radio constraint, the maximum number of traversed routing paths 
for any data source node is limited to be |iV|-2, where |iV|-2 represent the number of sensor nodes in the 
network besides the sink node and the data source node. In other words, at Step 7, the times of 
negative decision to go back to Step 3 is bounded to be \N\-2. Hence, we do not have an infinite 
looping problem. 

At step 1 of Figure 5, the link arc weight for multicast group g on each link / is setting to be: 

X u 2 gdl 

Arc si = dJ rU i + u l> +a ' + ^ u )ph + 2 Uj<t 9 + x X X u llk s jkPll . (18) 
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The first two terms are from the subproblem 3 where the physical meaning for these two multipliers 
is the tree constraint violation cost for multicast group g. Note that I do a summation of u 2 gdl with all 

the data source nodes of multicast group g. By this summation, we could ensure the link arc weight is 
the same for every data source node d. As a consequence, the union of the routing paths for every data 
source node becomes a tree. The third term is the power transmission cost. The forth term is cost of 
violating the transmitting and receiving channel assignment. 

The fifth term is the cost of violating the co-channel interference for the termination node of link /. 
If j" is the termination node of link /, then ^ ^ '^ u lk £ jkPij = X X M </'<t £ /'<t • ^ u ifk £ fk calculates the 

ieF jeNkeN ieFkeN keN 

channel violation cost of the sensor nodes that are within sensor node j" interference range 
(i.e., Sj-k = 1) that use the same channel i. By summation over all the available channels, ^^ ul„ k 8 n 

ieFkeN 

is the cost of violating the co-channel interference for the termination node of link /. In other words, if 
the termination node of link / violates the co-channel interference constraint, the ^ ^ ^ u yk £ jkPij 

ieF jeN keN 

cost will be high enough such that this link / will be unlikely to be selected again in the next iteration 
of the Lagrangean solution process. To summarize, by this link arc weight setting, I try to identify the 
multicast tree not only from minimizing transmission power point of view but also considering the 
channel assignment and co-channel interference. 



Figure 6. Illustrative example for LGR-Primal algorithm. 
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An illustrative example is given in Figure 6. In Figure 6, each sensor node is assumed to be 
equipped with two radios. The transmission cost on each link is assumed to be 1. In Figure 6(a), it 
shows the final data aggregation tree and the channel assignment when co = 3. If we have only two 
available non-overlapping channels, then data aggregation tree in Figure 6(a) is not a feasible solution. 
According to Step 5 in Figure 5, all the five incident links to a node G will be assigned with arc weight 
Z such that the new routing path will skip a node G. Similarly, the next routing path traverse via a node 
F or a node H will also be infeasible due to channel constraint. This is because for the node F or node 
H, there are two incident data source nodes such that there is no feasible channel for the node F or 
node H under two available channels. So all the incident links to the nodes F and H will be assigned 
with arc weight Z such that the next routing path will skip the node F or node H. Finally, we could 
have a feasible data aggregation tree that traverse a node / as indicated in Figure 6(b). Note that if a 
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node I is taken away from Figure 6, then we will not have a feasible solution since the total arc weight 
on the data aggregation tree will exceed Z as shown in Step 7 in Figure 5. In this case, after traversing 
three routing paths (i.e., traversing via G, F or H), we can conclude that there is no feasible solution 
when a node I is taken away from Figure 6. 

The computational complexity for the GIT algorithm is 0(\N\ ) for each data source node. The 
channel assignment (i.e., step 3 of Figure 5) for any sensor node on the data aggregation tree needs to 
identify the channel assignment for the other sensor nodes that are within its interference range. Hence, 
the computational complexity is 0(\F\\N\) for each data source node. Note that, for any data source 
node, the maximum number of iterations (Step 7 to Step 3 in Figure 5) to identify the channel/radio 
constrained routing paths is limited to be \N\-2, where \N\-2 represent the number of sensor nodes in the 
network besides the sink node and data source node. In other words, at Step 7, the times of negative 
decision to go back to Step 3 is bounded to be \N\-2. The computational complexity is 0(|F||7V| ) for 
each data source node. 

In the following, I show the complete algorithm (denoted as LGR) to solve Problem (P). The 



computational complexity of the above LGR algorithm is o(\F\ \og\f\ + ^|Z) 



for each iteration. 



Algorithm 1: LGR Algorithm. 

Begin 

Input: Network topology, data source nodes 
Output: Data aggregation tree 

Initialize Lagrangean multiplier vectors u' (0) = 0, Vz = 1,....,6. 

UB = |G|xX a / an d LB = a very large negative number (e.g.,-oo) //upper and lower bounds, 

/el 

respectively 

quiescence _age = 0, and step_size = 2. 

For iteration = 1 to Max_Iteration_Number, perform the following: 

Solve subproblem 1, subproblem 2, subproblem 3, subproblem 4, subproblem 5. 
Compute Z LR in (LR). 

If Z LR {u)>LB 

LB = Z LR (u) and quiescence _age = 0. 
Else quiescence _age = quiescence age + 1 . 
If quiescence _age = Quiescence Threshold 

step_size = step_size/2 and quiescence _age = 0. 
Run LGR-Primal algorithm (Figure 5). 
Compute the new upper bound ub. 
If ub < UB then UB = ub. 
Update the step_size. 
Update the Lagrangean multiplier vectors. 
End For 
End 
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6. Computational Experiments 



The sensor nodes are randomly placed in a 1 x 1 area. Sensor nodes are uniformly distributed on 
the deployment area. If there are 100 sensor nodes, the deployment area will be partitioned into 100 
grids with the same size and each sensor node is placed in the center of the grid. The most top left node 
is selected as the sink node such that we could have a data aggregation tree with larger depth. The 
transmission cost is equal to the square of the Euclidean distance of the link (i.e., signal attenuation 
constant a = 2). Two different types of data source nodes are simulated. The first one is event-driven 
where neighboring sensor nodes beside the event will become the data source nodes. The second one is 
random-source where data source nodes are determined in random. Hence, the data source nodes in 
event-driven are clustered but the data source nodes in random-source are randomly placed in the 
deployment area. 

For LGR, Max_Iteration_Number and Quiesceince Threshold are set to 1,000 and 30, respectively. 
The step _size _coefficient is initialized to be 2 and will be halved when the objective function value of 
the dual problem does not improve for iterations up to Quiesceince Threshold. The experimental 
results are within minutes of computation time. All of the experiments were running in a PC with 
INTEL™ P4 1.8 GHz CPU. In order to evaluate the solution quality of LGR algorithm, I compare 
LGR with the other four algorithms proposed in [19] under different varieties of parameters. 



Figure 7. Performance comparison with respect to traffic loads. 
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In Figure 7, it shows the results under increasing number of multicast groups with respect to loose 
channel and radio constraints. There are two important observations. First observation is that the LGR 
algorithm is superior to the other four algorithms with respect to total transmission cost. Second 
observation is that in random source model, the solution quality of LGR is constantly superior to 
ICADAR in random source model for both light and heavy traffic demands. In event-driven model, the 
LGR outperforms ICADAR especially in heavy traffic demands. Recall that in event-driven model, the 
data source nodes are clustered so that it will consume more radios and channels as compared to 
random-source model. In event-driven model, the channel and radio constraint will be more stringent 
at heavy traffic demands (e.g., 9 multicast groups). This indicates that to incorporate penalty cost for 
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violating channel assignment and co-channel interference constraint on the arc weight, the selected 
data aggregation trees are resilient to more stringent channel and radio constraints than the other four 
algorithms. 

Figure 8. Performance comparison with respect to the channel and radio in event-driven. 
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In Figure 8(a), it shows the transmission cost with respect to number of available channels under 
loose number of radios. It is expected that the total transmission cost will be increased at small 
available channels due to more stringent channel constraint. It is interesting to observe that, for SPT, 
GIT and CAGIT algorithms, the total transmission cost remains the same when <x> is decreasing. It is 
because that they belong to an one shot algorithm, which will not alter its routing decision with respect 
to different a> . On the other hand, ICADAR and LGR algorithms could choose its routing decision not 
to aggregate under stringent channel constraint. Such kind of behavior is also observed in Figure 8(b). 
In addition, it is interesting to observe that stringent radio constraint incurs even significant increasing 
in transmission cost as compared to stringent channel constraint. It is because the maximum number of 
predecessor (i.e., incoming links) for each sensor node is limited by the number of radios on the sensor 
node. Then in stringent radio constraint, it results in longer data aggregation tree which makes it 
difficult to perform data aggregation. LGR algorithm is superior to the other four heuristics in both 
Figure 8(a) and Figure 8(b), and it is significant in stringent channel and radio constraint. 

In Figure 9, it shows the performance comparison with respect to the channel and radio in random- 
source. As compared to Figure 8, the solution quality of LGR over the other four heuristics is more 
significant. In Figure 9(a), besides LGR could get much lower transmission cost data aggregation tree, 
LGR could locate feasible data aggregation tree under stringent channel constraint (4 channels). We 
also observe the same performance improvement of LGR in Figure 9(b). As compared to ICADAR 
could only locate feasible data aggregation tree with at least 4 radios, LGR could locate the feasible 
solution even with 2 radios. By carefully examine the difference between Figure 8 and Figure 9, we 
could conclude that LGR outperforms the other heuristics under stringent channel and radio constraint, 
especially in random-source model. 
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(b) Limit of Radio vs. Total Transmission Cost 



Institutively, in the dense WSN (i.e., large network size), the number of feasible data aggregation 
trees will also be increasing. In Figure 10, we could observe that the total communication cost is 
decreasing with respect to increasing network size. Recall that transmission cost is defined to be the 
square of the Euclidean distance of the link. From the objective function point of view, a concatenation 
of multiple links with shorter distance is preferable than a single link with longer distance. However, 
co-channel interference will be more significant such that it is not easy to identify feasible solutions. 



Figure 10. Performance comparison with respect to the network size. 




When network size is above a certain threshold, we could not get feasible solutions due to the co- 
channel interference constraint will not be satisfied. LGR and ICADAR algorithms could identify 
feasible solution even when the number of sensor nodes is 196. In addition, LGR could locate lowest 
cost aggregation tree. This reveals that with capturing the penalty cost of co-channel interference to 
construct the multicast tree, LGR algorithm is superior to the other four heuristics under variety of 
network size. From the comparison between event-driven and random-source model, it indicates that it 
is not easy to identify another routing path for a randomly distributed data source node under stringent 
channel constraint. As compared to co = 8 in Figure 10(a), we have co = 10 in Figure 10(b). In other 
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words, we need more channels to have feasible solutions for sparse and dense sensor network in 
random-source model. 

In Figure 11, it investigates the total transmission cost under different communication radius. It 
could be expected that with increasing the communication radius, it increases the probability of co- 
channel interference. Hence, it will be more difficult to identify the feasible channel assignment on the 
data aggregation tree. It might need to identify another routing path not to violate channel constraint 
(i.e., longer routing path that is out of the interference range of the other paths) under larger 
communication radius. We could observe that LGR algorithm outperforms the other heuristic 
algorithms in event-driven and random-source. As compared to ICADAR algorithm will rapidly 
increase the total transmission cost, we have more mild increase total transmission cost for LGR 
algorithm at large communication radius. This indicates that by incorporating penalty cost for violating 
the co-channel interference on the arc weight of each link, LGR could identify more cost efficient data 
aggregation tree to meet the co-channel interference under large communication radius. 



Figure 11. Performance comparison with respect to communication radius. 
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Table 1. Performance comparison between LGR and the other four heuristics. 



Heuristic 


Network load 
(No. of groups) 


Channel 
Constraint (co) 


Radio 
Constraint 


Network Size 


Communication 
Radius 


SPT 


(200%, 200%)* 


(175%, 150%) 


(400%, 350%) 


(300%, 142%) 


(180%, 175%) 


GIT 


(200%, 125%) 


(175%, 100%) 


(400%, 250%) 


(16%, 96%) 


(180%, 83%) 


CAGIT 


(50%, 125%) 


(125%, 100%) 


(300%, 250%) 


(16%, 142%) 


(100%, 175%) 


ICADCR 


(31%, 13%) 


(17%, 25%) 


(50%, 100%) 


(7%, 7%) 


(61.4%, 92%) 



*: (Event-driven, Random-source) 



From Figure 7 to Figure 11, we observe that LGR algorithm outperforms the other four heuristics. 
Besides more efficient transmission cost, LGR algorithm can locate feasible solutions at high network 
load, stringent channel/radio constraint, dense sensor network, and stringent co-channel interference. 
We define an improvement ratio = ((L-M)/M x 100%)to be the performance metric, where M stands 
for the total transmission cost from LGR and L stands for the transmission cost from the other four 
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heuristics. When there are no feasible solutions for the other heuristics (e.g., in Figure 8(a) no feasible 
solution when the number of multicast groups is larger than 6 for SPT, GIT and CAGIT), L stands for 
the largest feasible X-axis value from LGR and M stands for the largest feasible X-axis value from the 
other three heuristics. Note that in Figure 9(a), smaller co indicate stringent channel constraint, L stands 
for the smallest feasible X-axis value from four heuristics and M stands for the smallest feasible X-axis 
value from the LGR. I summarize the computational experiments in Table 1 . The first number in the 
parenthesis is LGR improvement ratio in event-driven and the second number is in random-source. 

7. Conclusions 

Data aggregation that could eliminate redundant data transmission is particularly useful in the 
limited power WSN. However, data aggregation also incurs collisions such that it produces extra 
energy loss from retransmission. By assigning different channels to sensor nodes that are within each 
other's interference range, it could eliminate the problem of retransmission. This requires sophisticated 
data aggregation routing and the channel assignment strategies. Besides channel assignment, sensor 
equipped with multi-radios could transmit/receive from multiple sensor nodes simultaneously to 
minimize the latency. However, channel and radio are limited resources in the WSN that need to be 
planned carefully to minimize the total transmission cost without violating the co-channel interference 
constraint. This paper studies the channel and radio constrained data aggregation routing problem in 
the WSN. I model this CRDAR problem as a mixed integer and linear programming problem and 
proposed Lagrangean relaxation technique (LGR) to tackle this problem. Unlike the existing heuristics 
(SPT, GIT, CAGIT and ICADAR), the data aggregation tree is based on the transmission power. The 
proposed optimization-based heuristics could identify the data aggregation tree from the perspective of 
transmission power and channel/radio resources simultaneously. The proposed LGR outperforms the 
other four heuristics in terms of total transmission cost with respect to all kinds of traffic load, 
available channels/radio resources constraints, network size, and communication radius. 
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